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The concept of new methodology of adding QCD NLO corrections in the 
initial state Monte Carlo parton shower (hard process part) is tested numerically 
using, as an example, the process of the heavy boson production at hadron-hadron 
colliders such as LHC. In spite of the use of a simplified model of the process, all 
presented numerical results prove convincingly that the basic concept of the new 
methodology works correctly in practice, that is in the numerical environment of 
the Monte Carlo parton shower event generator. The differences with the other 
well established methods, like MC@NLO and POWHEG, are briefly discussed 
and future refinements of the implementation of the new method are also outlined. 

PACS numbers: 12.38.-t, 12.38.Bx, 12.38.Cy 
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1. Introduction 

Successful operation of the Large Hadron Collider (LHC) at CERN is result- 
ing in rich harvest of experimental data. Even more data at higher energy and with 
higher statistics will be available over the next two decades from the LHC exper- 
iments. One of the challenges in the proper understanding and interpretation of 
these data, possibly leading to discovery of new phenomena, will be perfect mas- 
tering of the "trivial" effects due to multiple emission of soft and collinear gluons 
and quarks. Perturbative Quantum Chromodynamics (pQCD) HUH [3], together 
with the clever modelling of low energy nonperturbative effects, will be the basic 
and indispensable tool for disentangling the Standard Model physics component 
in the data. 
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of the European Union Grant PITN-GA-20 10-2645 64 (LHCPhenoNet), the U.S. Department 
of Energy under grant DE-FG02-04ER41299 and the Lightner-Sams Foundation. 
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2. Overview of the method 

Most of the methodology used in this work is described in ref. 0. In the fol- 
lowing, additional details relevant for Monte Carlo (MC) implementation and tests 
are described. Let us start with a description of the initial condition of the forward 
evolution (necessary in the MC implementation) which was omitted in ref. (4). For 
the proper understanding of this implementation it is necessary to recall some ba- 
sic facts about the use of maximum rapidity of emitted partons (angular ordering) 
as the evolution time variable. 

While the energy of the emitted gluon is a natural variable to handle infrared 
singularities, the angular variable is best suited for controlling collinear singular- 
ities. The logarithm of the angle of the emitted gluon (rapidity) with respect to 
the emitter partem emerging from the initial hadron, is a natural "master variable" 
for modelling collinear singularities. The angular variable is also well suited for 
modeling the structure of the non-abelian soft limit (colour coherence) 0161171. 

Conversely, a hard process in deep inelastic lepton-hadron scattering (DIS) or 
Drell-Yan (DY) process acting as a "probe", either backscattering (in the Breit 
frame for DIS) or absorbing (into a heavy boson in DY) the emitter parton, in 
a well defined rest frame of the hard process (RFHP), has its own energy scale 
used also as a master variable in collinear factorization and renormalization group 
equations. More precisely, for the hard process with the center-of-mass energy Q 
( y/S in the DY process), the parton entering hard process has energy equal Q/2. On 
the other hand, in the same RFHP, the initial hadron energy is E% and the Bjorken 
variable is the ratio x = Q/(2Eh) (x is invariant with respect to boosts along the 
emitter direction). The luminosity distribution of this parton D(Q,x) is commonly 
referred to as parton distribution function, PDF in short. It is weakly dependent on 
Q and is measured experimentally at each value of Q separately, that is at a given 
value of Q varying the energy £/, = Q/ (2x) of the initial hadron seen in RFHP. 

The important practical question for Monte Carlo modelling of the emission 
of the collinear gluons is: how to relate the variable t = ln(<2/A) governing the 
"evolution" of the PDF in the traditional DGLAP scheme^] such as MS scheme, 
and the rapidity variable of emitted gluons? 

In the following, to answer the above question we shall consider for simplicity 
the LO case with non-running as, and for pure gluonstrahlung (non-singlet QED- 
like component of PDF) from a single emitter^ For the MC purpose, we define 
the evolution variable ( as a hypervelocity of the Lorentz boost from the initial 
beam hadron rest frame to RFHP, t = S (or equivalently, the hypervelocity of the 
beam hadron in the RFHP). For each emitted gluon we define the rapidity in 
the rest frame of the initial hadron, or q in RFHP. Next, we require in the context 



Formulated typically in terms of exponentiation of the collinear singularities or using the renor- 
malization group equations, or both. 
2 One hemisphere in DY, or initial state cascade/ladder in DIS process. 
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of pQCD description of the gluonstrahlung that |; < in the RFHP (in the initial 
hadron frame ^ < S). In other words, t = S is a limiting rapidity for emitted gluons 
(or | = 0). Of course, e* = ^ = where m/j is hadron mass and s = 4E|. 

The role of perturbative QCD is to relate D(Q,x) measured in two experiments 
A and B with probes at the scales Qa and Qb, provided that Qa > Qb » mu- Two 
PDFs, D(Qa,x) and D(Qb,x), will differ because of the gluon emissions located 
in the additional phase space within the (Ea, S#) rapidity (angular) interval. Also, 
experiments A and B will use different RFHPs, connected by the Lorentz boost of 
the hypervelocity At = Sa — S# = In 

What is now the difference between the more traditional choice of the evolution 
time variable i = In ? of DGLAP and our preferred definition t = S = In ^ = 

A r mh | RFHP 

In (maximum rapidity of the emitted gluons)? When comparing two exper- 
iments with hard probes at the scales Qa and Qb, At = ^{Qa/Qb), while more 
phase-space conscious At = At — \r{xa/xb)- The "offset" ln(^A As) is formally of 
the NLO clas{J] and can be neglected within the LO approximation^] hence both 
choices are equally good at LO level. However, the use of (angular) t assures the 
completeness of the phase space of the emitted gluon, no gaps (nor "dead zones"), 
so it is the preferred choice in the MC modelling, aiming at the NLO level evolu- 
tion in the next steps. Additionally, the parallel use of i = t — In is quite useful 
and essential for other purposes, like introduction of the running as, etc. 

In particular, t is more natural for defining the initial point of the forward 
evolution (the stopping rule in the backward evolution). In order to assure the 
validity of pQCD it is required that the energy scale of the probe qo » A, nth is 
reasonably above the non-perturbative scales, like A ~ rtif, ~ lGeV, at the above 
initial point. This leads to the initial forward evolution point at to ~ ln(^o/A) and 
to ~ \n{qo/mh) — lnxo, as implemented in the following MC. In other words, gluons 
with rapidities below to are regarded as "unresolved", i.e. to = is a maximum 
rapidity for all unresolved gluons. 

It should be noted that the above discussion is quite standard in the context 
of any Monte Carlo parton shower using angular ordering. This line of the MC 
partem shower inspired by CCFM model [8], see also refs. |]9l [TO] [TTJ , is presently 
developed by CASCADE MC authors lfl2l . In particular, when using (maximum) 
rapidity t as the evolution time variable in the time ordered exponential of the 
QCD parton distributions the complete multigluon phase space is covered (with no 
gaps, "dead zones"), while the straightforward use of the ordering in the t variable 
in the MC would result in gaps between emitted real hard gluons, see also brief 
discussion of the corresponding kinematics in ref. |[T3l . 



3 It induces extra 0(0$) term in the evolution kernel. 

4 We have to remember to take it into account at the NLO level, when defining NLO evolution 
kernel. 
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2.1. Single LO ladder — basic building block in the MC 

Let us define multigluon distribution in the single initial state ladder taken 
in the LO approximation, which is a building block in our parton shower MC 
implementation, as an integrand in the following "exclusive/unintegrated PDF": 

D(t,x) = J dx dZ 5 x=Xo z d (t ,xo) G (t,t Q - In x Q \Z), 

G(t,to\Z) = e-^i(f\[ d*<E{h) 6^ (1) 

n=0\i=\ J 71 / 

where P(z) = \{\ + z 2 ), to = ln(^o/A). The "eikonal" phase space integration 
element is defined as in ref. \4r\ 



j3wi , d 3 k 1 dtydk+ je 
d 3l E (k) = — K -=■ =jt-i— d% 
y ' 2k Q k 2 271 k+ ^ 

and k^ = k° ±k 3 . In the above we use rapidity t, = \ In jpr \~ defined in the beam 

hadron rest frame Rh, while r] = | In | rfhp °^ re ^ - ® was defined in laboratory 

frame J^] They are simply related by i; = In ^ — r\. Rapidity ordering is now t = 

^max >£«>•••> & > >•••> £o = *o. where t = = ln(q /m h ) - lnx . 
The direction of the z axis in the RFHP is traditionally pointing out towards the 
hadron momentum. A lightcone variable of the emitted gluon is defined in the 

2k + 

usual way as a,- = of the emitter (after i emissions) as Xj =x$ — L}=o a j> an ^ 
finally z% = Xi/xi-i. The Sudakov formfactor Sf is determined by the "unitarity" 
condition 

l 

dZG(t,t \Z) = l (2) 



o 



and we omit its explicit definition, which involves the usual cutoff 1 — n < 6 reg- 
ularizing the IR singularity ^ = 4^-. The above feature is instrumental in the 
Markovian MC implementation, which provides D(t,x) for any value of t > to. 

The initial distribution do(qo,xo) can be related to experiment, to previous steps 
in the MC ladder, or to PDF in the standard MS system. Its precise definition is not 
essential for the following tests of implementation of the NLO corrections to the 
hard process, hence we will define it only numerically. We only notice that due to 
eq. Q the baryon number conservation sum rule 

/ dxD(t,x)= / dxo do{to,xo) 
Jo Jo 



5 A single ladder (parton shower from single emitter) is defined in the "tanget space" of momenta 
k, see ref. |4). 

6 Or, alternatively, in the overall center of the mass system. 
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is preserved. 

Finally, note the use in eq. (n| of the rescaled four-momenta kf 1 within the 
"tangent space", as defined in ref. 141. The mapping W — > k M can be defined^] once 
the ladders are connected together with the hard process, back in the common 
standard Lorentz invariant phase space, see ref. [4] and the following sections. 



2.2. Two-ladder LO multiparton distributions 

As a necessary introductory step to correcting the hard process to NLO level, 
let us start with defining and testing our simplified MC parton shower implement- 
ing the DY process with two ladders and the hard process, all three in the LO 
approximation j^] 



G = / dx 0F dx 0B d {to,xo F )d (i ,x B) £ 52 / dxp dxB 

J ni =Q n 2 =0-J 

•'-< r ln 1 \ (= 1 JL / 



(3) 



xe- S B I I [ld-Uk,)% n , \ ^(^/) 

i>n-. 2 \;=i 

xdz 2 (P- £ kj\q u q 2 ) -j£(sx F x B ,Q)WjZff > . 
j= l 

In the LO approximation we set W^c = 1- This weight will be defined/restored 
in the next section. In the above we use rapidity variable r|, defined in the overall 
center of the mass system (CMS). Rapidity £, of eq. ([!]) is translated into T|, dif- 
ferently in the forward part (F) of the phase space r\op > r|,- > S where we define 
5/ = In ^ - r\ u and in the backward (B) part E > rj, > rio B where = - In & + rj, 
should be used. The boundary between the two hemispheres E is for the moment 
set to be at S = 0, but in a more sophisticated versions of the MC will be cor- 
related with the position of the produced heavy boson (LO), or heavy boson and 
the hardest gluon (LO+NLO)0For the initial condition in the evolution we define 
T|of = ln(<7o/"i/i) — ln(^of) and = — \n(qo/mh) +ln(jcofl)- However, for the 
sake of simplicity we will set \n(qo/mi r ) = in the following. 

Phase space integration of eq. ^ for = 1 and using eq. ([I}, provides us 

7 This mapping preserves the rapidity variable. 

8 In this work we adopted notation of ref. |4|, in particular d%2(P',qi,qi) = 8' 4 '(/ > — q\ — 

, cPqi d i q 1 



v 



Gluon phase space is always fully covered (no gaps). 
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with the classical factorization formula 

Cq= dxp dxB Df(t,Xf) DB(t,XB) Gb(sxf%b)- (4) 
Jo 

In testing numerically the above formula, the convolutions Dp(t ,Xf) = (do®GF)(t ,xf) 
and Ds{t,XB) = (do (g) Gb) {t ,xb) are obtained from separate simple Markovian LO 
Monte Carlo exercises. As was stressed in ref. 0, the above LO formula rep- 
resents our LO MC without any approximations and can be tested with arbitrary 
numerical precision. Such a precise numerical test is demonstrated in the next 
section. 

In eq. Q distributions are expressed (similarly as in eq. ([TJ) in terms of the 
kf four-momenta in the tangent space. The mapping V 1 — > W is understood to 
be exactly the same as defined in ref. El, that is done simultaneously for both 
hemispheres, using ordering in the variable jr],- — S|. The details of this mapping 
do not influence the validity of eq. ([4]). 



2.3. Two LO ladders and NLO- corrected DY hard process 

Introduction of the NLO corrections to the hard process is done using a single 
"monolithic" weight on top of the LO distributions of eq. In the fol- 

lowing numerical exercises we will implement defined exactly as in ref. 0]. 
Let us recall this definition in a slightly more compact notation, for the sake of 
completeness: 



W, 



NLO 
MC 



l+A 5+ y+£ 



Pl(<?1^2,^) 



+ 1 



(5) 



with the NLO soft+virtual correction Ay + $ 
tion part: 

(l-$) 2 do Bt , 



P(z Fj ) da B (s, Q)/dQ £ B P{z Bj ) do B (S, B)/dQ. 

! and the real correc - 



C F CL S 
7C 



:i* 2 



da 



(l-o) 



d ° B I- a \ 



l + (l-q-P) 2 do B . 
9a >P 2 dK {s ' Q) 



l + (l-q-P)2 do B fi 



d&a 



(6) 



The above represents the exact ME of the quark-antiquark annihilation into a heavy 
vector boson process with additional single real gluon emission^] and subtraction 
of the LO component already included in the LO MC. The angle 6 in the sub- 
traction (LO) part of the Born distribution is typically defined in the rest frame of 



10 We employ here the particular compact representation of ref. 1 14] of this ME as a combination 
of the Born differential sections with the redefined scattering angle 6. 
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the heavy boson, where q\+q 2 = 0, as an angle between the decay lepton mo- 
mentum qi and the difference of momenta of the incoming quark and antiquark 
6 = Z(<7i,/?of — Pob) P] while two angles in the NLO exact ME are defined pre- 
cisely as Qf = Z(qi,—pog) and Q B = Z(eji,pop). (The implementation of NLO 
corrections in POWHEG scheme in ref. lfI31 uses the same form of the exact ME.) 
Note that in the above we only need directions of the vectors pop and po B , which 
are the same as of the hadron beams. The variable s = sxpx B = (qi + q 2 ) 2 is the ef- 
fective mass squared of the heavy vector boson. Finally, we specify the lightcone 
variables a ; and p 7 of the emitted gluon for j in the F and B parts of the phase 
space: 

a j = l-z Fj , $j = aje 2( ^- s \ for j e F, 

h = l-ZBj, aj = $je- 2 ^- z \ for j E B. 

The above relations are explained in ref. H as resulting directly from the kinemati- 
cal projection operators extracting the LO part from the exact matrix element. Note 
that variables oij and p^- in the above relations are defined in terms of which do 
not obey the overall 4-momentum conservation. The transformation W- — > and 
its inverse (where k^ do obey 4-momentum conservation) are defined explicitly in 
the above work. Slightly improved (LO level) kinematical mapping, better suited 
for the NLO-corrected hard process will be proposed at the end of section 4. 

The exact phase space integration of eq. ^ including of eq. ^ is again 

possible, see ref. [4] for details, providing a compact expression for the total cross 
section: 

l 

ai = J dxp dx B dzDp (t,x F ) D B (t,x B ) Cb(szxfx b ) {& z =i{l +A S+V ) + C 2r (z)}, 


(7) 

where C 2r (z) = Z ^ [- 1(1 -z)] was derived in ref. 0. 

3. Numerical results 

In the following we shall first check that the simple formula of eq. Q with 
two collinear PDFs agrees numerically with the parton shower MC of eq. Q with 
the LO hard process (W^° = 1). Once the above "LO benchmark calibration" is 
successful, we shall check numerically whether the NLO formula of eq. ([7J) agrees 
with the MC integration of eq. (|3|, switching on the NLO correcting weight 
of eq. (|5]>. In both MC exercises we expect deviations only up to statistical MC 
error, or other imperfections of the numerical implementations. 

Other similar choices of the angle in the Born distribution are also perfectly valid within the 
LO MC. 
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Fig. 1: In the upper plot the LO distribution of rfo = ^ln(xp /xb) from the CMC 
LO parton shower (purple) and from the strictly collinear formula (green) are 
shown. The lower plot shows the ratio of the two, the agreement of < 0.5% is 
obtained. 



3. 1. LO benchmark 

Figure [^represents a "calibration benchmark" for the overall normalization at 
the LO level. In fact, we show in Fig. [T] the properly normalized distribution of 
the variable r]^ = j\n.{xp /xb)- In the collinear limit this variable represents the 
rapidity of W boson. This variable will differ substantially from the true rapidity 
of the W boson in the presence of the gluon with hard transverse momentum, 
but this is an acceptable approximation in the current LO exercise. The W mass 
distribution with the sharp Breit-Wigner resonance lineshape is not very interesting 
and we do not show it here. 

In Fig. [T] one of the distributions is from the MC generation of the variables 
(xf,Xb) according to 2-dimensional integrand of eq. This is done using the 
general purpose MC program FOAM ifToll . However, in this MC we need the 
collinear PDF D(t,x) in the entire range of x and t as an input. This distribution 
has been obtained from a separate high statistics run (10 10 events) of a simple 
Markovian MC (MMC), recording the resulting D(t,x) in the 2-dimensional ta- 
ble (a finite grid). In fact this MMC run solves the LO DGLAP equation (for 
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gluonstrahlung LO kernel) using the MC method, similarly as in refs. flllEODQ 
From the look-up table recorded during the MMC run, a simple interpolation is 
employed to obtain D(t,x) for any values of t and x in the next step, that is in the 
2-dimensional integrand used by FOAM. 

Another distribution in Fig. [T] comes from the full scale MC generation (four- 
momenta conserving) according to eq. The MC run with 10 8 events was used. 
In the MC implementation we cannot use the Markovian method because of the 
narrow Breit-Wigner peak due to a heavy boson propagator. We could employ a 
backward evolution algorithm of ref. GUI , but instead we have opted to employ 
a variant of the constrained MC (CMC) technique of ref. |[T3l . In fact, we com- 
bine two CMC modules and FOAM into one MC generating gluon emission from 
the incoming quark and antiquark which annihilate into the W boson. The LO 
hard process ME of the W boson production is implemented^ FOAM is taking 
care of the generation of the variables xp ,xb,xfo,xbo and the sharp Breit-Wigner 
peak in s = sxpxg, then initial parameters for two CMC modules are set and the 
gluon four-momenta W- are generated. Once they are mapped into k^, following 
the prescription defined in ref. 0, the overall energy-momentum conservation is 
achieved. 

Figure [I] demonstrates a very good numerical agreement between do/dr^ 
from our full scale LO parton shower MC of eq. Q and the simple formula of 
eq. (|4]) in the strict collinear kinematics (just convolution of two PDFs and the 
Born cross section). The LO MC is working in the standard phase space, with 
the exact 4-momentum conservation and agrees with precision < 0.5% with the 
simple collinear formula of eq. Q. The visible numerical bias is most likely due 
to finite size of the grid used to parametrize PDFs from the MMC run. 



3.2. Numerical test ofNLO correction 

Having cross-checked very precisely the overall normalization of our LO MC, 
we are now ready to do a similar cross-check in case of the NLO-corrected hard 
process. 

Figure [2] represents a principal (technical) test and proof of concept of our 
new methodology for implementing the NLO corrections to the hard process in 
the parton shower MC. The NLO correction to the r|Jj, distribution is obtained 
on one hand within the full scale parton shower MC featuring the NLO-corrected 
hard process as in eqs. Q and Q, and, on another hand, with a simple collinear 
formula of eq. ^ in which two PDFs are convoluted with the analytical function 
C2r{z), the "coefficient function" for the hard process. In Fig. [2] we present the 

12 The use of the MC method is not mandatory here - we could solve it using finite step methods, 
as in ref. 1191 . 

13 In the presend MC exercise the average over angular distribution of the W boson decay products 
is taken. This is irrelevant for the conclusions of our study and this averaging can be undone 
rather easily. 
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Fig. 2: The pure (— ) NLO correction to the distribution of rfo = \ \&(xf/xb) in 
CMC LO parton shower in W boson production (purple). It agrees with the strictly 
collinear formula (green) to within < 1% of the NLO correction itself. 



NLO corrections obtained using both calculations. The LO component, cross- 
checked in the previous section, is present in the MC but not shown in this plot in 
order to increase the "resolution". In Fig. [2] we also include the ratio of the NLO 
corrections from the two sources. 

As seen in Fig. [2} the result of the parton shower MC with the NLO-corrected 
hard process and the result of the simple collinear formula of eq. ([7]) agree very 
well, within the statistical error. 

In Fig. [2] we see only the NLO corrections, but how big is the NLO correction 
with respect to the LO? We show this in Fig. [3] where both the LO and NLO 
components are compared, and the NLO/LO ratio is plotted as well. As we see, 
the NLO correction to the rapidity-like variable for the W boson is only about 1.5% 
of LO, and this is unusually small. 

For this particular hard process the NLO correction is negative, hence in both 
Figs. [2] and [5] it is multiplied by the factor (—1), in order to facilitate visualization 
of the results. 

At the technical level, the inclusion of the NLO correction in the parton shower 
MC is straightforward, we are just activating of eq. In fact MC is pro- 

viding the LO and NLO-corrected results in a single MC run with weighted events. 
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Fig. 3: The comparison of the LO (purple) and of the pure (— ) NLO corrections 
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Fig. 4: The distribution of the NLO weight of eq. 



The weight is well behaved, strongly peaked near W^r = 1 , positive, and 

without long-range tails. The distribution of this weight is shown in Fig. |4] In the 
MC implementing the collinear formula of eq. (|7]), we again use FOAM, but now 
the generation space is 3-dimensional due to the presence of additional variable z. 
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Note that in all numerical results shown so far we have put Ay+s = 0, as it is 
completely unimportant for the purpose of the presented numerical analysis. For 
our study we also used toy initial distributions do(qo,t), which were parametrize 
as follows: 

xD qep (x, lGeV) = 1xu{x) +xd(x) + —xs(x), 

2xu(x) =2.19x 1/2 (l-x) 3 , 
xd(x) = 1.23x 1/2 (l-x) 4 , 
xs(x) = 1.35x 02 (l-x) 7 . 

4. Discussion and comparison with other methods 

The new method of introducing the NLO corrections in the hard process pro- 
posed in ref. lU and tested in this work is clearly very different from the well es- 
tablished MC@ NLO EQ and POWHEG E2H23 methodologies. Ref. offers a 
limited discussion on these differences. Having at hand MC numerical implemen- 
tation we may elaborate on certain issues in more detail, in particular we are going 
to show numerical results illustrating differences with the POWHEG technique. 

At first sight, the most striking difference with the POWHEG and MC@NLO 
techniques are: 

• "Democratic" summation over all emitted gluons, without deciding explic- 
itly which gluon is the one involved in the NLO correction and which ones 
are merely "LO spectators" in the parton shower. 



The absence of (1/(1 — z))+ distributions in the real part of the NLO cor- 
rections (kinematics independence of the virtual+soft correction). 



In the following, we shall elaborate mainly on the first point, analyzing in a 
detail how of eq. ^ is distributed over the multigluon phase space. In order 

to make the discussion maximally transparent, let us consider a simplified weight 

w^=i+^wr, wr= (8) 

f^F P(z F j) da B (s,Q)/da 

which is limited to one ladder (one hemisphere). Moreover we put it on top of 
the MC modeling gluon emissions from single quarkj^jessentially the multigluon 
distribution of eq. ([TJ. 

We start by examining the inclusive distribution of gluons on the Sudakov 



logarithmic plane of rapidity £, and variable v = ln(l — z). This is shown in Fig. 5a 



14 We use the Markovian MC implementation, but optional use of the CMC would provide iden- 
tical results. We use quite a wide range of /, corresponding to = 7TeV. 
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(a) LO 1-gluon distr. p(t,v). Plateau at 2C | = 0.10 




(a) 



0.08. 
O.O67 
0.04 \ 
0.02i 



(b) 



Fig. 5: (a) The inclusive distribution of gluons on the log Sudakov plane of rapidity 
t = ^max and v = ln(l — z). (b) Contributions from all gluons weighted with the 
component weight Wj . 



The distribution looks as expected, and the flat plateau represents IR singularity 
ICf^d^Y 1 ^ (for constant a^) with the drop by factor 1/2 towards z = 0, due to 



l+z 2 



It ' 1— z 

factor in the LO kernel. 



In Fig. 5b we show contributions from all gluons weighted with the component 
weig ht -Wf 10 of eq. Q. (We insert a minus sign in order to facilitate visualiza- 
tion.) Here we see that the NLO contribution is concentrated in the area near the 
rapidity of the hard process t = ^ max , which has to be true for the genuine NLO 
contribution. On the other hand, the fact that the NLO correction dies out towards 
the IR limit z — > 1 is not guaranteed in the collinear factorization. It results from 
the conscious choice that our LO differential distributions reproduce the correct IR 
limit not only in LO but also in NLO and in the entire phase spacef^] 

Another important point is the completeness of the phase space near the (z = 0, 
t = ^ max ) phase space corner. Both POWHEG and MC@NLO use standard LO 
MCs which feature an empty "dead zone" in this region^ which is critical for the 
completeness of the NLO corrections in the hard process. They have to fill in this 
empty part of the phase space with MC events according to the correct LO+NLO 
distribution. Correcting for this deficiency of the standard LO MC requires non- 
trivial effort. In our case, the problem of the phase space incompleteness is ab- 
senj^Jand we simply re- weight the LO distribution (MC events) to the NLO level. 



Older version of the standard LO MCs do not always reproduce the correct soft gluon limit 
beyond the LO level. 

6 This is due to the use of the boost transformation in the standard LO MCs to get the overall 
four-momentum conservation. We avoid this transformation (problem) as we use the rescaling 
transformation only in kj — > kj in the LO MC. 

7 Of course, reconstructing the LO parton shower also requires non-trivial effort. 
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Finally, we also see that the NLO correction is very small, which might be a 
general feature of the new method. It is mainly due to the absence of the (1/(1 — 
z))+ terms in the NLO correction - this is a separate issue discussed in ref. H, see 
also a few remarks below. 



LO gluon K=I 



LO gluon K>1 



0.06 
0.04 

0.02^ 





Fig. 6: The inclusive distribution of gluons of Fig. 5a split into the hardest (in k T ) 
gluon (left) and the rest (right). 



Looking at Fig. [5] it is tempting to conclude that the dominant contribution to 
Y,jWj W may come from the gluon with the highest ln/:J ~ ^ + ln(l — zj), that 
is the closest to the hard process corner (z = 0, t = ^ m ax)- We may easily relabel 
gluons generated in the MC, — > Y.K> such that they are ordered in the variable 
Kk = + hi(l — zk), with K = 1 being the hardest one (k.k+\ < K/j). 

In Fig. [6] we show a split of the inclusive distribution of Fig. 5a into the K = 
1 component (hardest gluon in k T ) and the rest K > 1. As we see, the K = 1 
component saturates/reproduces the original complete distribution of Fig. [5a| over 
all the region where the NLO correction (Fig.[5b|) is non-negligible. This is exactly 
the observation on which POWHEG technique is built! Moreover, as noticed by 
the POWHEG authors, taking the K = 1 component is sufficient to reproduce the 
complete NLO correction (up to NNLO). 

Let us check numerically the above statements by means of comparing the 
NLO correction to the x = xoJljZj distribution from Y.j^f L ° an d from W^f. 
This comparison is shown in Fig. [7] As we see the K = 1 component saturates 
the entire sum very well, whereas the K = 2 component is quite small. A natural 
question is: why bother to keep the entire sum instead of taking only the K = 1 
contribution? In fact we can, which is valuable feature of our scheme. However, 
we stress that in the POWHEG scheme the K = 1 gluon is generated in the MC 
separately in the first step and other gluons are generated (by the LO parton shower 
MC) in the next step. This is fine and easy if the LO MC uses & r -ordering, while in 
case of the LO MC with angular/rapidity-ordering additional effort of generating 
the so called vetoed showers and truncated showers is needed in the POWHEG 
method. In our method, the angular ordering is used but the vetoed/truncated 
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showers are not needed, even if we replace the sum 

Zj^f by K = 1 compo- 
nent W^lo . is there any rationale for keeping the sum over gluons in NLO weight 
at all? There are two reasons for keeping it, at least optionally: (a) the valuable 
crosscheck of the NLO MC against the simple collinear formula of eq. ([7]) is exact 
only if we keep the sum, (b) it may turn out that keeping the sum reduces missing 
NNLO corrections. In our opinion one should keep both versions and check which 
one better fits the complete NNLO or better agrees with additional resummations 
beyond LO. 

As an additional illustration for the above discussion, in Fig. [8] we show the 
distribution of gluons ordered in rapidity, starting from the gluon with the maxi- 
mum rapidity, the closest to hard process. As we see the gluon distribution with 
the highest rapidity £, ~ ^ max (7=1) features a ridge extending towards the soft re- 
gion. It is important to notice that the width of this ridge goes to zero when £ — > 
in the IR cut-off (1 — z) < £■ Hence, sooner or later the gluon with the highest 
£, will not be able to reproduce/saturate the gluon distribution in the NLO corner 
close to hard process, and the NLO correction will be highly incomplete. In the 
& r -ordering this is not the case, of course. 

In the above, we have mainly discussed the differences of our MC with the 
POWHEG method in which, similarly to our case, the negative MC weight is not 
allowed. The MC @NLO method roughly corresponds to generating in a separate 
MC branch events according to non-positive NLO correcting distribution. In the 
same MC branch, events filling the empty phase space near the hard process corner 
are also addedF^I 

Both MC@NLO and POWHEG feature the (1/(1 - z))+ components in the 



This luckily reduces the number of events with the negative weight. 
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ALL: p (x,v} = t s ^-; Plateau at 2C I = 0. 10 



Angular ordering: J=l 




Angular ordering: J=2 





Angular ordering: J=7 ~<n> 




Fig. 8: The distribution of gluons ordered in rapidity, as in our basic LO MC. 



NLO corrections, which are the source of the practical complications there, while 
they are absent in our approach. In MC@NLO and POWHEG case these (1/(1 — 
z))+ corrections act effectively as the "in-flight" translation of the PDFs from the 
MS collinear factorization scheme (FS) to the FS used effectively in the MC, (see 
ref. H|). We propose to shift this translation beyond the MC, as a rather simple 
redefinition of PDFs which should be done "off line", from the point of view of 
MC. The above issue requires a dedicated study (in preparation), and is also closely 
related to the upgrade of the ladder part of the MC to the NLO level. 

Having in mind that the considered method is more general, and can be also 
applied for introducing the NLO corrections in the middle of the ladder I?24ll25l . it 
is an interesting question whether limiting the sum Y,j W^ LO to one (or two) terms 
would/could be used in order to upgrade the QCD evolution in the parton shower 
to the complete NLO level, which would open many new promising avenues in 
the development of the high quality QCD parton shower MCs for LHC and other 
colliders. This question will be addressed in the forthcoming study in ref. |[26l . 

The present work provides a numerical crosscheck of the ideas outlined in 
ref. 0. We briefly mention the most urgent future studies which will necessarily 
follow this work (some are already completed but unpublished). The two most 
important issues are: (i) adjusting the choice of S at NLO level, and (ii) select- 
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ing a better choice/definition of initial PDF. Also, adding missing graphs for the 
NLO corrections, that is graphs with gluon to quark transitions is needed. This 
should be simpler than the presented gluonstrahlung contributions due to lack of 
IR singularities. 

The present choice of the rapidity boundary S = is good at the LO level, 
and it also correctly reproduces the integrated NLO cross sention. However, at 
the exclusive level, the exact NLO distributions must be properly reproduced in 
the limit when all gluons but one are collinear (have small k T ), for instance, the 
rapidity difference between the heavy boson and the hardest (in k T ) gluon. For the 
above aim the best choice is to identify S with the rest system of the heavy boson 
and the hardest gluon T|*. This can be easily obtained by means of refining the 
mapping k, — > k, in such a way that it is used twice. For the first time with S = 0, 
then the rapidity rj* is determined and kj — > kj mapping is repeated with E = r\*. 
Obviously some gluons will be reclassified as belonging to another initial beam 
ladder^ The above solution was already tested and works correctly. 

Concerning further refinements on the initial PDF, this issue would be resolved 
automatically if MC was fitted to the experimental data or if the PDFs have been fit- 
ted within the MC scheme. If the initial PDF is to be taken from a standard library 
of PDFs in the MS scheme, then it will be necessary to correct it using the differ- 
ence of the counterterms of the MS and MC schemes (see eq. (44) in ||4]]). From 
the classic analysis ll27ll of NLO corrections to the DY process, it is known that 
this correction will be large and dominated by the term ~ ^ ln ((^) ( j n the 

region where quark distributiorj^jis strongly varying in x. Note that in POWHEG 
method the above correction is implemented in the MC by means of explicit gen- 
eration of the variables in the convolution implementing NLO corrections and the 
corresponding manipulation on four-momenta is done. In contrast, in our method 
NLO corrections are included entirely through MC weight and no extra kinematics 
transformations are needed (beyond these of the LO MC modelling). 

In the above context an interesting numerical results are presented in ref. ||28l 
- they illustrate size and location of the x-variation in PDFs due to kinematics 
manipulations in POWHEG driven by NLO corrections. In our method the entire 
kinematical modification of the longitudinal parton fraction x is due to the LO 
mapping and shape modification due to the NLO weight. We expect this effect to 
be less sizable in our method, but a separate study would be needed to verify it. 



Luckily, this "flow" of gluons from one to another hemisphere does not influence the overall 
MC weight. 

Similar phenomenon will occurs for the gluon distribution. 
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5. Summary and outlook 

A new methodology of adding the QCD NLO corrections to the hard process in 
the initial state Monte Carlo parton shower is tested numerically using heavy boson 
production at hadron-hadron colliders. The ladder parts of the parton shower are 
modelled in the LO approximation, also using these new methods. The presented 
numerical results prove that the basic concept of the new methodology works cor- 
rectly in the numerical environment of a Monte Carlo parton shower. The differ- 
ences with the well established methods of MC@NLO and POWHEG are briefly 
discussed. Also, possible refinements of the method are indicated. 

Clearly the "proof of concept" is successful, and more work is required before 
practical application will emerge. 
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